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Abstract 

We study the radiation of gravitational waves by self-gravitating binary systems in the low- 
energy limit of Hofava gravity. We find that the predictions for the energy-loss formula of 
General Relativity are modified already for Newtonian sources: the quadrupole contribution 
is altered, in part due to the different speed of propagation of the tensor modes; furthermore, 
there is a monopole contribution stemming from an extra scalar degree of freedom. A dipole 
contribution only appears at higher post-Newtonian order. We use these findings to constrain the 
low-energy action of Hofava gravity by comparing them with the radiation damping observed 
for binary pulsars. Even if this comparison is not completely appropriate - since compact 
objects cannot be described within the post-Newtonian approximation - it represents an order 
of magnitude estimate. In the limit where the post-Newtonian metric coincides with that of 
General Relativity, our energy-loss formula provides the strongest constraints for Hoi'ava gravity 
at low-energies. 
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1 Introduction 

General Relativity (GR) continues to stubbornly agree with every observation related to gravity 
[1]. This would be extremely desirable if the theory could be merged with quantum mechanics in a 
straightforward way. Unfortunately, the current situation is far from this: the search for a consistent 
theory of quantum gravity remains elusive and there is no experimental guidance to shed light on it. 
Furthermore, the cosmological constant problem asidj^l, the success of GR as a low-energy effective 
field theory (EFT) points towards the Planck mass Mp ~ 10^^ GeV as the physical frontier where 
one expects to learn anything about quantum gravity. If the preceding arguments are realized in 
Nature, experimental information about quantum gravity will indeed be sparse in the foreseeable 
future. 
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'^One may argue that the cosmological constant problem is a hint towards the actual theory of quantum gravity, 
and that a successful framework of quantum gravity should provide a mechanism to explain this phenomenological 
observation. We do not address this particular issue here. 
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More interesting for phenomenology are the proposals for ultraviolet (UV) completions of GR 
where the previous logic fails. These include models of gravitation with a low-energy cutoff beyond 
which GR ceases to be valid If this cutoff scale is as low as the TeV, these proposals may have 

interesting phenomenology and may even be relevant for the resolution of the hierarchy problem. 
Another recent proposal in this category is Hofava gravity [H [5]. This framework yields a concrete 
UV completion of GR, with effects that may permeate basically any gravitational experiment. It is 
on the implications of Hofava gravity for gravitational radiation that we pursue in this paper. 

Essentially, Hofava's proposal consists of considering the existence of a preferred time-foliation 
of spacetime. Assuming the presence of this absolute structure, the GR Lagrangian can be supple- 
mented by operators which render it power-counting renormalizable without destroying the unitarity 
of the theory The result is a non-relativistic theory of quantum gravity [5] (in the sense that it 
is Lorentz-violating). The preferred foliation is in principle detectable at any energy scale, and it is 
not surprising that this approach (which is designed to cure the unsatisfactory behaviour of GR at 
distances of the order Mp^ ~ 10~^^ cm) generically also modifies the theory at large distancej^ [S]. 
Among the different implementations of Hofava's idea, we consider the so-called "healthy extension" 
[7]. This version possesses a stable Minkowski background where the issues about strong coupling 
appearing in other approaches are absent. Furthermore, variants of Hofava's original proposal can 
be retrieved for a particular limit of this (generic and stable) case [?]■ 

The low-energy (large-distance) sector of the theory is encoded into a scalar field ip, calleci§ 
the "khronon", that characterizes the foliation structure and interacts with a metric field. We 
refer to this low-energy scalar-tensor theory as "khronometric" theory 0. The extra scalar field 
ip turns out to be massless, and its presence modifies most of the predictions of GR, including 
the parametrized-post-Newtonian (PPN) parameters [5l[71[9] and cosmological phenomena |10| 
These modifications differ from those of standard scalar-tensor theories [H [12]. They are close, 
however, to the predictions of Einstein-aether theory (or ae-theory for short) |13| . This is not 
surprising since both theories incorporate a field whose expectation value violates Lorentz invariance 
(a unit timelike vector in the case of Einstein-aether, and ip in our case), and are otherwise generic. 
It can be shown that the khronon ip corresponds to the hypersurface-orthogonal mode of ^-theory, 
and many of the predictions of both theories are indeed identical [5l |8]. The PPN parameters 
derived from ^-theory and khronometric theory restrict the parameter space of those theories but 
are otherwise in agreement with current observations. Thus, both (low-energy) models represent 
interesting alternatives to GR, which, furthermore, have a high energy cutoff . The further advantage 
of khronometric theory is that beyond this energy cutoff there is a known UV completion, in the 
form of Hofava gravity. 

counterexample to this argument can be found in [6]. However, it is not clear how GR is recovered at large 
distances in this proposal. 

^From Greek xpofoc, - time. The khronon is also known as the "T-field" [8]. 
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The aim of this paper is to further constrain khronometric theory based on the loss of energy 
due to the emission of gravitational waves (GWs) from a binary self-gravitating system. This is 
a relevant test for gravitational theories given its sensitivity to the way gravity propagates (e.g. 
the degrees of freedom and corresponding properties), and also to the strong- field regime since 
known astrophysical sources of GWs tend to have strong gravitational self-energies |12l HH ES] • The 
confirmation of GR's famous quadrupole formula in the damping of a binary pulsar's orbit is indeed 
one of its triumphs |16| I17| . Radiation tests have also been used in the past to constrain possible 
modifications of GR [H [181 [IHI- A priori for both ae-theory and khronometric theory, one expects 
this radiation formula to be modified due to a different speed of propagation of the tensor modes 
and the presence of new propagating fields. These modifications imply new ways to constrain the 
parameter space of the theory, independently of PPN and cosmological considerations. While the 
above expectations have been verified for ^-theory in the weak-field regime in |20| . the constraints 
obtained are not final since the astrophysical systems for which radiation damping has been observed 
are not in the weak-field regime [1]. The incorporation of strong-field effects in the Einstein-aether 
began in |21| . 

We focus on the radiation formula in the post-Newtonian (weak field, slow-motion and weakly 
stressed |14j ) regime of khronometric theory. This restriction is interesting for two reasons. First, 
we find deviations from GR's quadrupole formula already at leading order. (This is similar to what 
happens in ae-theory, as computed in |20|.) Second, and from a purely pragmatic point of view, many 
of the formulae we present in this paper are useful for the phenomenologically relevant situation of 
compact sources. First results relevant for the study of gravitational radiation from these systems 
include the black hole solutions derived in |22| 123) , and those for neutron stars in ae-theory |24) . 
The use of binary pulsar observations to constraint Hof ava gravity was suggested in |25) . 

To extract information about the damping of the orbit of a binary self-gravitating system from 
the emission of GWs, we take advantage of the fact that khronometric theory is semi-conservative (in 
the language of |12|). Then, for the bound system there exists a conserved energy that decreases due 
to the emission of gravitational radiation. By computing the energy fiux at infinity, we can derive 
the fiux of energy lost by the binary. Under the assumption that this energy is extracted entirely 
from the orbital motion of the binary, the subsequent damping of the orbits can be computed using 
Kepler's third law. This assumption has been tested to lowest order in GR |15| . and is plausible for 
khronometric theory. 

This work is structured as follows. In Sec. [2l we define the action for khronometric theory and 
the equations of motion relevant for low-energy phenomenology. Sec. [3] is devoted to the linearized 
equations for the fields far away from the source (far-zone). In Sec. [U we study the conserved 
properties of the source relevant for the post-Newtonian (PN) calculation. We derive the explicit 
expressions for the the different waveforms, up to and including the first PN order corrections in 
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Sec. O In Sec. [UJ we determine the formula for the average power loss in GWs. This formula 
is evaluated for a Newtonian system of two point-masses in Sec. [71 where the Peters-Mathews 
parameters for khronometric theory are derived. We summarize our results and conclude in Sec. [51 
Appendix [XI contains a derivation of the PPN parameters for our model (whose full expressions 
appear here for the first time). Appendix [Bl compares the monopole contribution, or lack thereof, 
in both khronometric theory and ae-theory for a particular choice of parameters. Finally, Appendix 
[Q provides a summary of the notion of energy relevant for our study. 

Conventions 

We use the (H ) signature. For an arbitrary expression X, the overbar X denotes the part of X 

linear in perturbations. The superscript X^^ is the non-linear part of X, i.e. = X — X. The 

dot X denotes the derivative of X with respect to time. Greek indices refer to spacetime, whereas 
Latin indices refer to space only. Repeated Latin indices are to be summed, e.g. Xu = 5^^Xij. We 
define the symmetrization of indexes as T(^ij) = ^ {Tij + Tji). We choose units where c = h = 1. 

2 Action for khronometric theory at low energies 

As outlined in the introduction, Hof ava gravity is based on the existence of an absolute time foliation 
of spacetime. This allows for the GR Lagrangian to be supplemented with higher dimensional oper- 
ators that render the theory power-counting renormalizable [1]. These operators are suppressed by a 
scale whose magnitude is constrained by various phenomenological tests. The most stringent of 
these tests comes from absence of deviations from Newton's law at short distances [5] which implies 
that M* > (10 /xm)~^ ~ 10^^ Hz [Il[5]. Thus, these higher dimensional operators are irrelevant for 
the binary systems of interesi0 and we neglect them in the following. The presence of a preferred 
foliation also has consequences at energy scales below M*. Indeed, at low-energies new operators 
appear (compared to GR) that are compatible with the group of gauge invariance preserving the 
preferred foliation, i.e. the foliation-preserving diffeomorphism [H [5]. Renormalization group ar- 
guments imply that these relevant operators should be added to the GR action, which has been 
done in the Stiickelberg (or covariant) formulation of the theory in [SllM]- In this formulation, the 
preferred-time foliation corresponds to the expectation value of a scalar field ip called the "khronon". 
This field is such that the normal to the surfaces of constant field is timelike, 

d^^d^'ipX). (1) 

^As an example, the famous PSR 1913+16 binary pulsar has a characteristic frequency of 10^ Hz [T]. We assume 
that the speed of propagation of all the modes is similar to the speed of light. We comment on this assumption when 
we derive the energy-loss formula in Sec. [Sj 
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The action of the theory is invariant under diffeomorphisms, and Lorentz invariance is broken by 
condition ([T]) in a spontaneous way. Also, the action must be endowed with invariance under field 
reparametrizations 

V^^/M, (2) 

which follows from our requirement of a preferred time-foliation as opposed to a preferred time. It 
corresponds to the time reparametrization invariance of the theory in the original formulation of [3] . 
The invariance under the transformations ([2]) is readily implemented by making the action depend 
on Lp through the combination 

r). in 

(3) 



^JdpipdPif 

Clearly, is non-singular whenever condition ([T]) is satisfied. Notice also that is a unit timelike 
vector. 

The low-energy action for the healthy extension of Hofava gravity corresponds to the most 
general action describing the coupling of (p with a metric field g^jj at low-energies and compatible 
with the aforementioned invariances |5]. It is given by 

' (4) 



where M;, is an arbitrary mass parameter to be related to the Planck mass. 



and a, /3 and A are free dimensionless constantqj. We also introduce a term Sm in Eq. (j4]) rep- 
resenting the matter component of the theory. We assume that matter is universally coupled to 
the metric g^^^ which enforces the weak equivalence principle [T]. This action defines what we call 
"khronometric theory". For later convenience, we introduce 

M ^ /" /Vf ^ /" 



■'X — 2 / ^ V crp 
where 

is used to compactify notation. 

Khronometric theory can be considered on its own as an alternative to GR with an extra 
scalar field, independently of quantum gravity motivations. This approach is similar to the way 
Einstein-aether theories are constructed. The only difference is that the vector is taken to be a 
generic timelike vector in ^-theory |13| . meaning that it has more degrees of freedom than in the 
khronometric case. It also implies an extra term in the generic action with respect to Eq. (j4]). 



^Note that the parameter A corresponds to A' in the notations of [5]. It also differs from the A defined in [4]. 
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This extra term can be absorbed by the ones in action dU for hypersurface orthogonal vectors, i.e. 
whenever satisfies Eq. ([3]) . Khronometric theory and ag-theory share the nice feature of having 
a high energy cutoff. The advantage of the former is that a UV completion in the form of Hofava 
gravity is known. 

Let the khronon and matter energy-momentum tensors be, respectively. 

The explicit expression for T^^y reads 

+ a a^ay + IVpK^i^^Uy-j - Uf^Uyu'^V pK^^ - 2a aaUf^pV y^u" + a a^apUpUy, 
where we have introduced the notation 

ap = u^VpUp. 

The equations of motion derived from varying the action with respect to the metric are 

Qm. = Gpy - M-2 (T^x^ + T-) = 0. (5) 
The equation of motion for the khronon field is 

= VpJ^' = Vp (^-^V^" [^aK% - « V.uH j = 0, (6) 

where 

As usual, this equation follows from the covariant conservation of the khronon energy-momentum 
tensor. That it can be represented as the conservation of a current is a consequence of the invariance 
of the theory under reparametrizations of the khronon given by Eq. ([2]) |26] . 



3 Equations of motion in the far-zone 

The physical system of interest for radiation damping consists of an isolated self-gravitating astro- 
physical source. By this we mean that there is a region of spacetime surrounding the source where 
the fields acquire their background values plus small perturbations. Thus, there exists a coordinate 
frame where the metric in this region satisfied, 

®In this section, Greek indices are manipulated with the Minkowski metric. 
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with \hfj_j^\ <C 1. For the khronon field, we fix the parametrization invariance ([2D by working with a 
time coordinate corresponding to the background of the field. Thus, we expand it as 

where |x| ^ ^- It is easy to verify that the background fields are indeed solutions of the equations 
of motion ([5]) in the absence of matter. To derive the flux of energy lost by this system, it is enough 
to understand the behavior of the fields produced by the isolated source in this region where they 
are weak. This is so because the energy carried by GWs is radiated away and eventually permeates 
the "weak-field" zone. We can extract the power radiated by integrating the flux of energy over a 
sphere surrounding the source at a particular time after the emission. This calculation is further 
simplified in the region far away from the source due to the applicability of the both the "weak-field" 
and "far-zone" approximations (see below). 

To derive the equations governing the perturbations h^i, and x, we split Eq. ([5]) and Eq. ([6]) into 
linear and non-linear parts as follows 

G^, - = M- V, 2x = -2x''- (8) 

The expression for t^i, reads 

T,. = T;ru + {T^u) - M^G^,^ . (9) 

This separation into linear and non-linear parts allows us to solve for h^i, and x perturbatively in 
Mj"^. The terms T^jy and can be interpreted to be source terms for the linear equations at 

different orders in M^~^. They include contributions from both matter and non-linear gravitational 
fields of lower order. For this paper, we are interested in matter sources that are weakly self- 
gravitating, slowly movinj^ and weakly stressed. These are known as post-Newtonian (PN) sources 
[14| . For these systems, one has 



i;,V2| 
V ~ |/io[, I 







rpm 

ij 


rpm 
-^00 




rpm 
^00 



1/2 

« 1, (10) 



where v is the typical velocity of the source. Thus, we can introduce v as a new parameter of 
expansion and consider the predictions of the theory at different orders in v, also known as PN 
orders. We content ourselves with the first PN corrections, which amounts to considering Eqs. ([8]) 
where the source terms also include PN corrections. In particular, the metric should be substituted 
by its first PN expression (Appendix |A]) whenever it appears in the non-linear source terms. 

This straightforward analysis is only suited for the first PN corrections. Beyond that, the 
analysis becomes more complicated due to the presence of tails and retardation effects. The correct 



^For theories with modes propagating at different speeds, this means that the typicai veiocity v of the source is 
small with respect to all of them. 
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treatment of the problem in general involves the separation into a near-zone and a wave-zone. In 
the near-zone, one can find the metric to any PN order including non-linearities and minimizing 
retardation effects. This corresponds to an expansion in the small parameter to desired order in 
V. In the wave zone, one can solve the equations of motion perturbatively in the fields and match 
the solution to the one found in the near-zone in a region where both approximations are valid 
|14l [T^ [271 12S]- For the first PN corrections considered in this paper, this analysis reduces to the 
one outlined in the previous paragraph. For higher order corrections the matching is much less 
trivial [H [El EH [29] . 

The linearized khronon energy-momentum tensor satisfies the following conservations laws 



(11) 



It follows from the invariance of the linearized theory under linear diffeomorphisms. Next, from the 
transverse properties of the G^i, and when the equations of motion are imposed, one finds 

{TJZ + T^u - MiG^J^) = 0. 
Together with Eq. (|lip . this yields the following conservation equations for the source tensor r^,y 

(12) 

which is of particular importance in Sec. [4] and beyond. 



3.1 Wave equations 

We decompose the gravitational perturbations into irreducible representations of SO{3), 



hi 



2(j), h 



di 



Oi 



1 -0, 



where tu = ditij = diVi = djFj = 0. We also define the Laplacian hy A = didi. 



(13) 



Tensors and vectors 

To single out the tensorial part of the equations of motion as written in Eq. 
transverse-traceless projector Pij^kr and the transverse projector Pij 



F'ij,kr 



PikPjr 



— Pin Pi 



ij ^kr 1 



Pi. 



A ' 



A straightforward calculation yields 



Pij,kr Q, 



kr 



-rPi 



ij,kr 



we introduce the 



(14) 
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leading to the wave equation for the tensor modes 

{ct^dl - A)ty = -2M-^Pij^ksTks, (15) 

with Cf = 1/(1-/3) representing the speed of propagation of the tensor polarizations. This coincides 
with the results of ae-theory |30) . 

Consider now the vectorial part of the equations. Contrary to ae-theory |20| [30] this sector does 
not contain any propagating polarizations. Indeed, one finds 



(1-/3) 



PijQoj = ^"' A [Vi - F^) - M^^PijTjo = 0, 

PikdjQkj = ^^^A (v - - F,) - M~^Pikdjnj = 0. 



(16) 



The first equation represents a constraint and its time derivative yields the second equation (which 
follows from gauge invariance). For definiteness, we choose to work in the gauge 

Fi = 0, (17) 

which completely fixes the gauge freedom in the vector sector. 

Scalars 

The scalar sector of the theory is different from GR. In particular, it includes an extra degree of 
freedom. We choose to work in the gauge 

X = B = 0, (18) 

which completely fixes the gauge in the scalar sector. The choices (jl7p and (jlSp are referred to as 
the "unitary gauge". In this gauge, the non-redundant equations of motion derived from ([5]) and ([6]) 
are 

aAcf) = 2Ail) - M-Voo, (19a) 
(/? + A)A^ = -2(A + 1)AV' + M^'^diToi, (19b) 

The speed of propagation of the scalar perturbation is given by 

2 ^ (a-2)(/3 + A) 
a(/3-l)(2 + /3 + 3A)' 

which coincides with the scalar mode of ae-theory 
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3.2 Far-zone expressions and post-Newtonian approximation 

The equations of motion (jlSp . (jl6p and (jl9p contain two types of equations that we wish to solve, 
Poisson and wave equations. The Poisson equation is of the form 

A^{t,x) = -47rp{t,x), 

whose solution for vanishing boundary conditions at infinity is given by 

j3~ P(*>*) 



\x — x\ 

We assume that the isolated source of GWs can be confined within a sphere of radius R. At distances 
far away from the source, r = \x\ ^ R, we can perform the expansion (f* = x^/r) 

\x - x\ =r -f'-x' + rO{R/r)'^. (21) 

We refer to this zone as the "far-zone". The leading contribution of the solution to the Poisson 
equation at large distances is then 

1 /"^s. 



^f{t,x) = -J d-'xp{t,x). (22) 

The sourced wave equations are of the form 

{c-''di - A) ait, x) = 47r/x(t, x), (23) 

with speed of propagation Ca- The solution to this equation with radiation boundary conditions is 
given by (see, e.g. |15) ) 

r,^ Kt-\x-x\/c.,x) ^ 

/ \x — x\ 



Besides adopting the far-zone approximation and using (|2ip . we also assume that r is such that 
r S> (jjR? /Cfj, where oj is the largest characteristic frequency of the source. This allows us to write 
the leading contribution as 

af{t,x) = - / d^x/i(t-r/c<^-Ff*xVc<x,x) = 5o / x p{t - r / , x) [f' x'/ c„)'^ , (25) 
r J r n. j 

where the last identity holds formally. This expression can be simplified further for the PN sources 
of interest \12\ [T^ I31| . As seen in the previous section, khronometric theory involves two speeds 
of propagation, the tensor and the scalar speeds q and c^. We assume that the system is slowly 
moving with respect to both speeds which are considered to be of the same order, q ~ ~ 1. 
Thus, for a typical velocity v ~ ujR of the source, the sum in Eq. ([25|) represents a well-defined 
expansion in the small parameter. 



i.e. it is a PN expansion, cf. Eq. (jlOp . In other words, every time derivative in the near zone 
represents an extra 0{v). 
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4 The source: conservation properties 



The source terms for the equations (|15|) , (|16|) and (|19|) are expressed in terms of the pseudo-tensor 
T^jy. In order to find solutions to the Poisson and wave equations in the far-zone, Eqs. (|'25p and (|2'2p 
indicate that we need to evaluate various integrals of Tfj_y. In what follows, we present results that 
are relevant for simplifying those integrals (and therefore the wave forms that appear in Sec. E]) and 
include leading PN corrections. We refer the reader to Appendix |X] for more details on the first 
order PN approximation and the parametrized-post-Newtonian (PPN) formalism |12) . 
From Eq. (|12|) . one can establish the useful integral conservation laws, 

j d^xTij = d^xfoox^x^' -^j d^xdf'ff.ox'x^ (26a) 
d^x foiX^ = - j d?x Tij . (26b) 
d^xfoox' = - / d^xTio+ / d^xd^Tf,ox\ (26c) 



In deriving the previous equations, we assume that all the boundary integrals cancel (the correc- 
tions to this assumption are negligible at large r). The difference with respect to the GR integral con- 
servation laws is the presence of the terms proportional to d'^T^o coming from the non-conservation 
of Tfj^i,, Eq. (jl2p . Remember that the current r^o is conserved for khronometric theory . Naively one 
expects d^T^Q to contribute to order as low as 0{v^). To see that this is not the case, we notice that 
Eqs. (j26p can be simplified by writing the equation of motion ([6]) as an equation for a conserved 
current (which corresponds to the Noether current related to the invariance of the theory under 
reparametrizations of Eq. ([2])), 

d^.{^gJ^)=Q. (27) 

Furthermore, in the unitary gauge, if = t and = (see Appendix D of [5]). Since J* is linear in 
perturbations, we find that 

Qx = -driV^JY"-- (28) 

Thus, 

d^^To, = Mid.i^jY"^ ~ 0{v'), (29) 

and the dipolar corrections turn out to be large in PN order. In particular, at order 0{v^) only the 
Eq. (|26cp is modified with respect to GR. A straightforward but tedious calculation using the PN 
metric displayed in Eq. (j48p yields 

I d'xd>^To^x' = \l'^'^P [(«™ - ^2'''')yr'' + ar'^W,™] + 0{v<'), (30) 
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where 

PP^ _ 4{a - 2(3) 
^ /3 - 1 ' 

PPN _ {a- 2/3)(-/3[3 + /3 + 3A] - A + a[l + /3 + 2A]) ^ ' 

" (a-2)(/3-l)(/5 + A) 

These constants are the PPN parameters related to the violation of Lorentz invariance of the theory 
(see Appendix [A]) . In the limit of small parameters they coincide with those found in [5] . The 
potentials Vf^^ and VK^^^ are also defined in Appendix [Aj Finally, the form of the Eq. (j30p is 
identical to the one found for ae-theory |21) . 

The previous formulae (j26p and (|30p can also be derived by relating the pseudo-tensor r^j^ to a 
conserved (but asymmetric) object. Indeed, from Eqs. p2|) and ([6]) it is evident that 



satisfies 



V = 0- 



This object has a contribution linear in the fields. To build a quadratic conserved pseudo-tensor it 
is enough to add the conserved current found in (j27p and consider the object 



The resulting integral conservation laws for this object are then identical to the ones found in 
\12\ [2T] . The existence of this conserved quadratic current is a generic consequence of the theory 
being semi-conservative in the language of |12) . This conserved current satisfies 



Then, one can use the Eq. (4.103) in [12j to compute (|29p . Even if this method may save a lot of 
computations, it is inconvenient since the result in |12) is derived in the PPN gauge, whereas we 
are interested in the result in the unitary gauge ([30|) . 



5 Wave forms in the far-zone 

We are now ready to compute the explicit form of the wave solutions in the far-zone, which we do 
consistently up to 0{v^) in the PN approximation. For the tensor and vector modes, the solutions 
of Eqs. dUD, (HID, (I25D and ([26]) are (in the gauge Fi = 0) 

1 - •■ 1 - • a 

tij{t,x) = - „ Pij,ksQks{t - r/ct) - 7—772 Pij,ksr''Sks,a{t - r/ct) + 0{v ), (32a) 

47rM^r lirM^ct r 

r2 /If 



V^{t,x) = ( - / d'i P.jT,o{t,x) ) , (32b) 
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where 



Q{t)ij = I{t)ij - ^6ijlkk{t), Iij{t) = j dxToo{t,x)x''x\ 
Sks,ait)= / d^x Tks{t,x)x'^. 



The quantity Qij represents the quadrupole of Too- Note that in the far-zone, the longitudinal 
projector Pij of Eq. (I14p can be substituted by the longitudinal part of the algebraic projector, 



This substitution is valid up to 0{R/r) terms. The object Pij^ks is defined as 

Pij,kr = PikPjr '^PijPkr- 

Anticipating the results of Sec. [H we notice that the energy- loss formula depends on the time 
derivative of the fields. For the vector part, the previous expression yields 

= = (7 / "^'^ ^^J-^iO^*'^)) • (33) 

From the conservation law (fT2|) . this term can be expressed as the integration over the boundary 
of the transverse component of the source, which cancels away from the source, and we can neglect 
the vector perturbations altogether. 

Concerning the scalar field, from the wave equation (|19cp one finds 

+ — 1 d^x foo(t - r/cs,x)x'- + ^^^f'-f^f^ 1 d^xToo(i - r j Cs,x)x!-x^ x^ (34) 



where 



+ -r^SkkAt - r/cs) - ^^^j§^^r^r^r>^S,,At - r/c,) | + 0{v'), 

Cs Cs{P + \) J 

_ (/3-l)(ar^-2ar^) 
^ = 3(a - 2/3) ■ (^^^ 

Notice that the conservation law foo = 5j(roi + Ji) has been used to show that the first moment 
of Too is constant in time and therefore ignored in (j34p . From the results in the previous section, 
we see that the modification to the GR results appear at order O(f^). Notice also that it follows 
from ([30|) and the constancy of the integral of Tjo that the dipolar contribution in (|34p is 0{v^) and 
suppressed by the PPN parameters. Finally, the quadrupole terms in the tensor and scalar sectors 
differ slightly, as they depend on different retarded times. For the remaining scalar fields cj) and 
from Eqs. (jl9p one finds 

<t^ = -^ + TT7T— I roo, E = -^^iilv-. (36) 
a AirMfar J p + X 
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6 Energy- loss formulae for post- Newtonian systems 



The definition of the energy carried by gravity waves is non-trivial (see [32j for a review on the 
concepts of energy and momentum in GR). For the problem at hand, we follow the procedure of 
|21j (see also |33) ) and use the notion of energy for asymptotically flat spacetimes derived in |34j . 
Given an isolated source, we can compute the time variation of this notion of energy by performing 
an integral of the flux in the far-zone, which we idealize as being infinitely far away from the source. 
We associate this energy loss to the energy carried away by gravitational radiation. As shown in 
|34tl35|. this alternative approach is equivalent to the one based on pseudo-tensors used in standard 
computations of energy loss due to gravitational radiation [12^ [T5| I31| . We give a brief review on 
this method in Appendix ICl 

In deriving the energy-loss formula, we make the following assumptions. We start by assuming 
that our system consists of an asymptotically Minkowski spacetime at early times, with the following 
fall-off properties in the unitary gauge, 

g^.u = ll^.u + 0{l/r), d^g^, = 0{l/r''), x = 0. (37) 

As for the matter fields, we assume that they vanish asymptotically to ensure that there are no 
boundary integral contributions. The previous conditions allow us to define a convenient notion of 



conserved energy £, Eq. (j6ip . as the conserved charge associated to the invariance of the asymptotic 
solution under asymptotic time translationj^ . To compute the fiux of gravitational radiation, we 
consider the moment of time when the emitted GWs are already at spatial infinity, which means 
that the fall-off properties of the fields change to 

/i^~0(l/r), /i^~a;0(l/r), ~ a;/c,0(l/r). (38) 

The quantity £■ with these boundary conditions is in general divergent. Nevertheless, its change due 
to the radiation emitted during a finite interval of time is well-defined |36) . We focus on computing 
the time variation of Eq. (|64p . As shown in Appendix [Cl S- is finite for conditions (|38p and only 
has contributions that are quadratic in the fields. 

We also consider the time average of the quantity £ over several periods of the source, {£). 
The final averaged energy-loss formula is a standard observable in GW experiments (including the 
binary system of interest, where the observed damping of the orbits occurs after several periods) 
and the final expression is simplified since total time derivatives vanish when integrated. 

The final expression is further simplified after one takes into account the following considerations. 
From the form of the solution of the tensor modes, Eq. (j32ap . and of the field tjj, Eq. ()34p . then in 



^Even if this symmetry is broken by the background for the field ip, it is still a symmetry due to the reparametriza- 
tion invariance of the theory ((2]). 
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the far-zone these fields satisfy the equation 

C(j dia = — f'^a, (39) 

for the corresponding speeds of propagation. Remember also that in the far-zone, the tensor modes 
tij are transverse with respect to the algebraic projector r^tij = 0. For the vector part, we already 
showed that it does not contribute to {£), as its time derivative cancels, cf. Eq. ([33]). Similarly, 
the fields E and (j) always appear under a time or a space derivative. Thus, we notice that the 
dependence on the source appearing in Eq. ()36p will be either higher order in R/r for the space 
derivatives (and therefore negligible), or of the form 

J d^ifoo = J d^x = Ml j d?xd,{^jY^ = 0. 

So, only the -0 contribution for E and (j) is non-zero, and therefore these fields satisfy relation (j39p . 
In fact, the latter relation is also satisfied by the scalar part of hij. 

The previous considerations (and some algebra presented in the Appendix [C]) yield the final 
result for the rate of energy loss of the system, 

{£) = -^i d^r'(-k,U,-^-^^i^^). (40) 
4/5^ \ct ' ' acs I 

Whereas the radiation emitted in the tensor modes always decreases the energy of the system, the 
behaviour of the emitted scalar modes depends on the parameter a. We see that the emitted energy 
is positive for the range < a < 2, as expected since these values are also required for the stability 
of the Minkowski background (absence of ghosts) [?]• 

Up to this point, we have consistently worked to first PN order (which corresponds to and 
includes 0(t)^) in the wave-forms). Given the time derivatives in Eq. (j40p . substitution of the 
waveforms (j32ap and (pi]) yields the energy loss of the system from gravitational radiation up to 
and including 0(t)^^), although corrections already appear at leading order, 0{v^^\ To simplify 
what follows, we therefore focus on Newtonian sources and the corrections at this order. Substituting 
the waveforms to lowest PN order in the previous expressions and performing the angular integrals, 
we find the energy-loss formulaic 



where (recall Eq. ([35 



ct 2c,(a-2)' 4c,(a-2)' 



'in this final formula, we compute the quadrupole and monopole terms at a time when radiation from both 
the tensor and scalar modes simultaneously reaches the boundary of the isolated system. For different speeds of 
propagation Ct and Cs, the discrepancy in emission time is irrelevant for stationary production of GWs. 
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The final expression, Eq. (|41|) . differs from tlie GR result in two ways: the coefficient corresponding 
to the quadrupole depends on the parameters of the model; there is a monopole contribution already 
at this first Newtonian order. This is similar to the formula derived for ffi-theorj0 |21| . Let us note 
something quite remarkable in the context of the khronometric theory that we are studying. All of 
the Solar System tests are passed in the limit |af^^| ^ 1, |a2|^^^ ^ which can be achieved by 
the single requirement \a — 2/3\ <^ 1, cf. (|31|) . In this limit, Z = 1, the dipole term (l3Up cancels and 
the monopole contribution in Eq. (|41|) is still present. This last result contrasts with the ^-theory 
case for which there is only a modified quadrupole (in the equivalent limit). This discontinuity 
between the two theories is discussed in Appendix iBl 



7 Energy loss by a Newtonian binary system 



To complete the calculation, the power-loss formula ()4ip must be supplemented by the equations of 
motion of the system to desired post-Newtonian (PN) order. We content ourselves with a 2-body 
Newtonian system composed of point-masses mi and 771,2. The matter action is then given by 



2 

E 

A=l 



niAdsA , (42) 



where dsA represents the proper-time of the A-th particle. A priori, niA depends on the khronon 
field. Since we are only interested in the Newtonian system, it is enough to Taylor expand the mass 
around its background value and use only the leading order contribution, hence tjia is taken to be 
constant. Using the preferred time as the affine parameter, the energy-momentum tensor derived 
from (02]) is 



where the A-Xh body follows the trajectory x\{i) with four-velocity li^. At Newtonian order this 
yields 

2 

'iij{t - r/ca) = <9o ^ mAx\it - r/ca)x'^{t -r/c^). 
A=i 

We evaluate this at very late times as explained in the previous section. Next, from the geodesic 
equation derived from (|42p . we find Newton's law 

-i _ ^"^2.4 ..i _ ^ f^i 

Xi — — Gjv— 2~''l2' ^2 — Lt7V^~'^12' 
^12 "^12 



^"^For the monopole and quadrupole contributions, the definition of Z in [21] diflters from Eq. (|35p by a factor of 
two. We attribute this difference to a typo in the final formula for ip in |21) . 
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where we introduce (also see Appendix 

""^^ A.Mil-ay -12 = -!--- (43) 

As usual in binary systems, it is convenient to define the problem in terms of the relative distances 
and the position of the center of mass xqai = "^mi+m2^^ ' Finally, assuminj^ that the system is at 
rest with respect to the preferred frame (so that xcm = 0) we get 



r 



12 

with fi = mim2/M, M = mi + m2 and = is related to the expansion parameter v. Thus, the 
loss of energy in gravitational radiation for a Newtonian binary system is given by 

from which we deduce the 'Peters-Mathews' (PM) parameters |121 137| {hd = 0), 

Ki = 12(1 - a/2)^, = (1 -q/2) ^11^- 

Once the energy loss for the binary system is known, one can use Kepler's third law to relate it 
to the damping of the orbit. The expression for the change of the orbit's period for generic PM 
parameters in terms of other orbital parameters of the system can be found in |38j . 

In GR, the previous analysis suffices to predict the radiation damping of binary systems for 
compact (relativistic) sources, like the PSR1913+16 |12| . This is because the structure of the 
compact stars of the binary does not influence the orbit in GR (this is called the 'effacing principle' 
which is a consequence of the strong equivalence principle). This is certainly not true for most 
alternative theories of gravity. Thus, to yield concrete predictions about the radiation damping of 
systems with highly relativistic sources (sources with large self-energies), one must first understand 
the behaviour of the fields beyond the PN approximation. One can then use Eq. (j4ip to derive 
the energy loss resulting in a change of the orbit (at corresponding PN order). For scalar-tensor 
theories, the final result is a test of the strong-field regime [39 1 140 | HT]. For ^-theory, the first steps 
were performed in |21) based on the stellar solutions of |24| and using the effective field theory 
methods of gS] SSI El SS] 

It is beyond the scope of our article to derive the radiation damping of these realistic systems 
(including relativistic self-gravitating objects) for khronometric theory. In any case, we do not expect 
the new corrections to cancel the ones we have already derived for the Newtonian source, and thus 
we find it appropriate to use Eq. (j44p to set order of magnitude bounds on the free parameters of the 



Corrections to this assumption are considered as higher order in the PN expansion. 
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khronometric action Q. Current data on the radiation damping of the Hulse- Taylor binary system 
agrees with GR up to a level slightly better than one part in one hundred [U [T7] . This means that 
the formula (|44p should agree with GR to O(10~^), which finally implies the bound (for the case 
where a, /3 and A are of the same order) 

a~/3~A<10-^ (45) 

The previous bound is less stringent than the bounds coming from the PPN analysis [5l [12] , 

i«r^i < lO"^ ia™i < 10-^. 

As can be directly seen from Eq. (|3ip . the PPN bounds are automatically satisfied in the limit 
a = 2(3. In this limit, Z = 1, and our expression (|45p yields the most stringent bound for the 
theory. Notice in particular that it constrains the propagation speeds to be close to c = 1. Another 
constraint in this limit comes from the difference between Gn as derived in (|^3|) and the value for 
Newton's constant appearing in Friedmann's equation, Gc [Z]- The value of Gc is constrained by 



nucleosynthesis and satisfies 



Gc ^ 



< 0.13 [46], which, in terms of the parameters in the action 



([4]), implies the estimate a, f3, X < 0.1 [7j. Also, we should ensure the absence of gravitational 
Cherenkov radiation, which implies cf > 1 and > 1 (this means that a particle moving through 
the aether does not radiat a. Notice that the speeds are superluminal, which does not pose a 
threat to Lorentz violating theories as long as causality is maintained. 



8 Discussion 



Our aim has been to study the radiation loss from an isolated source in the PN approximation for 
khronometric theory. This theory is an interesting alternative to GR with a high energy cutoff and 
for which a UV completion is known in the form of Hof ava gravity. It is also very similar to ^-theory, 
as in both cases there is a preferred time coordinate. The difference is that khronometric theory 
has only one extra scalar degree of freedom, the khronon, whereas ^-theory relies on a timelike unit 
dynamical vector leading to three extra degrees of freedom, consisting of one scalar and one vector 
field. 

For arbitrary parameters, we have shown in Eq. (j4ip that the formula controlling the power 
loss of the system (which may be related to the change of the orbital period of a binary source) 
is modified with respect to GR already at lowest, Newtonian, order. In particular, the quadrupole 
contribution differs from GR, partly due to the different speeds of propagation of the tensor modes 
in both theories. Furthermore, there is also an extra monopole contribution at this order. The 

^■^We thank D. Levkov and S. Sibiryakov for pointing this out to us. 
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monopole at leading order in khronometric theory contrasts with the usual situation in other scalar- 
tensor theories |18| . At higher order, there are other modifications, including the dipole term (|3Up . 
Quite remarkably, in the phenomenologically interesting limit where all PPN parameters coincide 
with GR (which corresponds to the limit a = 2(3 for our parameters), the monopole is still present, 
and its strength is proportional to the parameters appearing in the action of the theory, Eq. 
These results for khronometric theory are similar to those of ae-theory, modulo vector propagating 
degrees of freedom that are absent for the khronometric case. There is a key difference, however, 
since ^-theory only has a modified quadrupole to lowest order in the equivalent limit. 

This work has been devoted to PN sources. These types of sources do not correspond to the ones 
found in the binary systems of interest, which are compact and characterized by strong gravitational 
fields. Despite this, we have evaluated the energy-loss formula for the simplest possible system: a 
Newtonian binary. Doing so provides an order of magnitude estimate on the parameters of the 
theory (as we do not expect corrections due to strong-fields to cancel the modifications apparent in 
the power loss formula). Thus, our results are relevant for constraining the case a = 2/3. In this 
case, requiring the rate of radiation damping to be close to GR sets constraints on this parameters 
of order O(10~^). These constraints represent the strongest phenomenological bounds for this 
particular choice of parameters and are relevant for the cosmological implications of the theory, 
including the recently suggested model of dark-energy |47) . 

Sources with strong self-energies is left for future research and can be treated in our theory in 
the same way as scalar-tensor theories |39| HOl I41| : a phenomenon of "scalarization" modifies the 
orbit of these sources as compared with the post-Newtonian ones. It would also be interesting to 
consider our results in the parametrized post-Einsteinian framework introduced in |48] (see also 
|49| for the binary pulsar constraints for this framework). Finally, the consequences of alternatives 
theories of gravity for experiments of direct detection of GWs have been recently discussed, see e.g. 
[H [50]. We hope to extend these works to khronometric theories in the future. 
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A Post-Newtonian expressions 



The PPN formalism is a valuable tool for comparing theories of gravitation with each other and 
with experiment in the weak, non-relativistic hmit [12] . In this section, we briefly present the steps 
involved in the PPN calculation for khronometric theory (see also [5]). The final result are the 
parameters (all the other PPN parameters cancel) 
pPPN ^ ^PPN ^ 1^ 

PPN 4(a-2/3) 

«i = ' (46) 

PPN _ (g - 2/3)(-/3(3 + /3 + 3A) - A + a[l + /? + 2A]) 
"2 - (a_2)(/3-l)(/3 + A) 

Notice that the PPN parameters for khronometric theory for arbitrary values of the parameters in 
(HI) appear here for the first time. They coincide with results in [5] in the limit of small parameters. 
The non-zero parameters af^^ and af"^^ indicate that khronometric theory violates Lorentz in- 
variance. These same two parameters are non-vanishing for ae-theory, although the dependence on 
the parameters a, /3 and A is different. In both theories, however, the relationship between af^^ 
and a2^^ is the same 



„PPN 
"2 



af^^ (2/3-a)(3A + /3 + a) 



2 (A + /3)(2-a) • 

To compute the previous results we closely follow \12\ [51] to which we refer the reader for further 
details. The source is assumed to be a fluid with a covariantly conserved energy-momentum tensor 



T^^ = (p + pU + p)v''v'' - pg''", 

where is the four velocity of the source, p the rest mass energy density, 11 the internal energy 
density and p the isotropic pressure of the fluid. The source is assumed to satisfy (llOp . 
In what follows, recall that the different fields have the following expansion, 

goo = l + 0{v^) + 0{v^), go, = 0{v^), 
g,j = -h, + 0(t;2), X = 0{v'') + 0{v^). 



(47) 



Also, we use the following potentials 

J \x-y\ 

where Gn is defined in Eq. (|43|) and the correspondence F ^ f \s given by 



U^l, ^i^ViVi, ^2^U, ^>3^n, ^i^p/p, 



yppN ^ ]yppn ^ '"ji^j ~ yj)(x' - y*) 

The steps to take are: 



\x — 
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1. Solve QoQ to order 0{v'^). For this we use the 00 component of Eq. ^ to 0(t)^), which yield 



2 

2. Solve gij to 0(7;^). Following |51| . we choose the gauge conditions 

The arbitrary constant T will be chosen to write the result in the PPN gauge. Then from the 
ij component of Eq. ([5]) to 0{v'^), we find 

2 

Ahij = 87rG]ypSij. 

3. Solve X to 0{v^). The khronon equation of motion ^ to leading order is given by 

Ax - Tdohoo j = 2(A + /3) 

4. Solve (7oj to 0(^3). In our gauge, the Oi component of Eq. ([5]) to 0(t;3) yields 

. 3 87rGNpv^{a - 2) - [-2 + a + r(l - PWihoo 
Aho. = j—^ . 

5. Solve goQ to 0{v^). From the 00 component of Eq. ([5|) to 0{v^), we find 

4 2 2 2 2 

A/loo = dihoodihoo - /iooA/ioo - 4A$i + 4A$2 - 2A$3 - 6A^>4 



+ 2/3(3 + /3 - 2r) + 2(3 + 3/3 - 2r)A + 2a(/3(r - 1) + (F - 3)A)) ^ ^ 
+ (a-2)(/3 + A) 



where H = —Gm f d^yp\x — y\ is known as the superpotential. 

6. To go to the PPN gauge, we choose F that cancels the term depending on H in the previous 
equation |12| . 



Putting everything together, we have (to desired order) 

gPP^ = l- 2U + 2U'^ - 4$i - 4$2 - 2^3 - 6$4 = 1 + AH, 
gf/"" = -(1 + 2U)6,j = 6i, (-1 + AH) , 

^PPN I ^PPN ^PPN\^rPPN I I ^PPN\ijrPPN 

%i =2^" + "! ~"2 Wi +2^1 + "2 Wi 

PPN ^ {a-2p){2 + ^ + 3X)H 
^ ~ 2(a-2)(/3 + A) ' 



^■^We use a number over the field to keep track of tlie order in v. 
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which, compared to the generic PPN metric (see for example, Eq. (A. 11) of |51| ) imphes that all 
the PPN parameters vanish except for the ones cited in (|46|) . 

The PN metric in the unitary gauge of Eqs. (|17|) and (|18|) is easily derived from these PPN 
expressions. It suffices to go from the PPN gauge to the unitary gauge via a diffeomorphism 
6x^ = S,^ satisfying 

^ ~ 2(a-2)(/3 + A) ' ^ ~ 2(/? + A) 
This leads to the following PN metric in the unitary gauge 

goo = l + AH + 0{v^), 



9,, = (-1 + AH) - i^±A±l^d,d,H + Oiv^), 

P + A 



(48) 



X = Oiv*). 

B The Einstein-aether and the monopole 

In both khronometric and Einstein- aether theories, we compare the monopole contribution to the 
energy-loss formula in the limit for which the PPN parameters are identical to GR. The free pa- 
rameters of khronometric theory are a, B and A and those of the Einstein- aether |13) are Cj for 
i = 1, . . . , 4. We have the correspondencq^ ci = 0, C2 = A, C3 = /3 and C4 = a. Notice that one less 
parameter is needed to define khronometric theory. This is because the action of a hypersurface 
orthogonal aether (which is equivalent to khronometric theory [5]) contains a term that can be 
absorbed by the others, reducing the number of independent terms from four down to three. 

Comparing the results of this paper and the work presented in |20) , we see that the waveforms for 
the spin-0 and spin-2 modes are essentially identical. The main difference comes from the expression 
for Z of Eq. (j35p (remind also footnote [T0|) . Let Z be the equivalent expression in ^-theory, 

3(C14 - 2ci3) 

where af^^, a^^^ are the Lorentz violating PPN parameters in ^-theory and Cij = Ci + cj. Then 
the khronometric expression for Z is precisely Z, but with ci = 0. 

The limit af^^ = af^^ = in khronometric theory can be achieved by setting a = 2/3 
and leads to Z = 1. By inspection of the energy- loss formula (j^Tj) . we see that the monopole is 



^*Recall that we are using the mostly minus signature. The mostly plus signature is used in |20) and leads to a 
different correspondence between the parameters. 
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proportional to Z and therefore persists in this hmit. In generic ^-theory, the equivalent limit that 
sets the PPN parameters to GR is given by 

C2 = — C4 = 50 

3ci Cl 

and leads to Z = 0. The corresponding monopole is proportional to Z and subsequently vanishes 
in this limit. Therefore, the values of Z and Z explain the presence or absence of the monopole in 
the limit when the PPN parameters are identical to those of GR. 

It is natural to ask if Z can be tweaked so that ae-theory has a monopole when of = a^^^ = 
0. A first possibility would be to consider the limit that ressembles khronometric theory, namely 
Cl = and C4 = 2c3. This leads to Z = 1, like in khronometric theory, indicating that a monopole 
may be possible. However, requiring only ci = implies that af^'^ = 8. One could try to set ci = 
and C3 = to get af^^ = 0, but this case of ^-theory has yet to be studied |13| . Alternatively, 
one may try to make the denominator in (j49p vanish to retrieve a finite limit. Setting C14 = 2ci3 
yields Z = 1. However, the second condition in ([50|) implies ci = C3 = C4 = 0, which is a singular 
limit for ae-theory. 



C Notion of energy for an asymptotically flat spacetime 

To characterize the energy carried away from a system by GWs, we use a method different from the 
standard technique defined in terms of the Landau-Lifshitz or related pseudotensors |12| [T5| EH [32] . 
Here, instead of computing the energy carried by GWs, we derive the loss of energy of the isolated 
system during the process of gravitational radiation. This resembles the definition of energy loss 
by the time variation of the Bondi-Sachs mass |281 132] . However, we will use a different notion of 
conserved energy that, to our knowledge, was first used in the context of GWs in |20| . This energy 
is well-defined for asymptotically flat spacetimes satisfying the boundary conditions (j37p . which we 
use to define isolated sources. Its conservation follows from the invariance of the asymptotic solution 
under time translations and it reduces to the standard notion of energy for flat spacetime |34| (see 
also |36( I52|). Since the method is not standard, this Appendix is devoted to presenting a succinct 
summary. We encourage the reader to consult the original literature to complement it. 

Given a Lagrangian density L{^) depending on some dynamical fields $, we define its associated 
4-form (we present the 3 -|- 1 case) as 

L(^>) = L{^)d^x. 

After integration by parts, the first variation of the previous form following from the variation 6^ 
can be expressed as. 
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where E$ = are the equations of motion of the theory. If the variation 5$ is a diffeomorphism 
generated by a vector field ^, the previous variation should correspond to the action of this trans- 
formation over L($), 

where L refers to the contraction of the form L with the vector field ^. Define the Noether current 
3-form associated to ^ and L{^) as 

Jl = 0l($,<55$)-%L. (51) 

This form is clearly closed when the equations of motion are satisfied. In practice, to find the 
components of the 3-form 0^, notice that it is dual to a 1-form. In components 

where the index of is risen with the metric g^'^ and e^^j^p are the components of the Levi-Civita 
3-form defined for the metric g^y. From this definition it follows that 

de^ = V^V^e^d^x = d^{^Ql)d^x, (52) 

from which one can easily identify the components of 0^ . 

To associate the flux generated by ^ to a Hamiltonian evolution from an initial hypersurface 
S, one must assume [341 152) that in the boundary of the initial hypersurface, denoted by 9S, it is 
possible to find a 3-form Bj, such that 

JdT. JdT. 

If such a current exists, the flux generated by ^ corresponds to the orbits generated by of the 
Hamiltonian 

H^= [ Jl- [ hBl. (53) 

Finally, since Jl is closed when the equations of motion are satisfied, it follows that locally Jl = 
dQL- Thus, when the equations of motion hold, can be written as a pure boundary term, 

Hi= [ (Ql-hBl). (54) 

To define a canonical notion of energy, we shall now assume that ^ is an asymptotic time 
translation, with components .^'^ — t- 6q and that the asymptotic conditions on the dynamical fields 
have been specified in such a way that the surface integrals appearing in Eq. (|54p approach a finite 
limit. The Hamiltonian then corresponds to the generator of time evolution. We define the canonical 
energy at a hypersurface slice of constant time Sf to be |34) 

£l = j^^ (Ql - Bl) , (55) 
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where represents the boundary sphere at the boundary of S^. Whenever £i is weU-defined, it is 
a conserved quantity, and we can remove the t label in Sf. 

We now apply the previous formalism to our action ([!]). The hypersurface of constant time 
corresponds to a sheet of the preferred foliation. Even if not necessary, it is convenient to work with 
an action for which 

/ k@L = ^- (56) 

This equation is not satisfied for the Einstein-Hilbert action part of ^ (see e.g. Eq. (87) in 
|34|). As explained in |34| . the existence of a background metric r]^y makes it possible to build a 
covariant action (which is required to get a conserved current (f5T]) ) equivalent to Einstein-Hilbert 
and satisfying (j56p . Indeed, let us write g^y = r]^y + h^y and consider h^y and rj^y as independent 
dynamical fields. We can then add a boundary term invariant under diffeomorphisms to the action 
(gD to yield 

s'^s + ^ I d^x(V^((r^,-fp5^^-(r;:,-r;:j5^")),,= J d'xL', (57) 

where T^u refers to the connection compatible with the background metric i]^y. The part corre- 
sponding to GR reads 

A/fS r 

The equations of motion derived from varying the previous action with respect to h^^y and r]f^y are 
the same, as these fields appear only in the combination Qf^y, except in the boundary term. As 
a consequence, rj^y can be considered to be Minkowski, and we can assume that the equations of 
motion fix h^y. 

For the computation of Jrr corresponding to the action (j58p and the vector field dt, with 
components 6q, we first notice that dt is a Killing vector of 7]^y, 

and the boundary term in Eq. (|58p does not contribute to Jrr- For the first term one finds the 
corresponding current 




This term is linear in the connection and does not depend on the derivative of 5g^y. To construct 
the conserved current, we use 
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Thus, under the assumption that the fields fall-off at large distances as (|37|) . the current (|59p vanishes 
asymptotically as 0(r~^), which means that its contribution to (|56|) cancels. Indeed the cancellation 
of the contribution to (|56|) holds in the more general situation where one considers variations 5gf^i, 
which do not change the asymptotic behaviour (|37p . Finally, the energy £'pr derived from Eq. (|58p 
coincides with the ADM mass which also agrees with the energy derived from the Landau-Lifshitz 
pseudotensor |lj4| [55] . 

The term in the action (jU yields a current 



{aa^V^u'^ - VpKP^) —5x + K^'-^d^^bx 



(60) 



X 



Remember that the invariance under diffeomorphisms is non-linearly realizeq^ on x 

kx = C° + edpX- 

From Eq. ^ this means that Jax ~ 0(1). Similarly = -^f^ ~ 5^0 + 0{l/r). Thus, 9: 
0{r~^), which means that the contribution of this term to ()56p cancels. 

Finally, we find that the conserved energy (|53p for the action (|57p inside a constant time hyper- 
surface Sf is given by 

£ = [ d^x^jl,, (61) 
with J'S, representing the coordinates of the 1-form dual to the corresponding 3-form, Eq. (j5ip . 



j^, ^ (ej<r + - 6^,L'. (62) 

The contribution from the khronon action is simplified once one considers the equation of motion 
for X- Indeed, in Eq. (j60p includes a term 

{aa^V^u'^ - VpKP^) ^ = r, (63) 
V A 



where J'^ is defined in ([6]). In the unitary gauge, this current is purely spatial, which means that 
this term does not contribute to (I6ip . 

We are eventually interested in the flux of energy loss through GWs, so we want to compute the 
quantity, 

S= j d^x^jl, = -(f dn^r^r'J's,, (64) 



'We could also work with the field tp for which S^ip = ,^''9,i<^. 
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where we have used the fact that the current Jg, is conserved on-sheh, which is a consequence of J 
being closed and (|52|) . The final ingredient is to evaluate Jg,. From Eq. 



rr ——^'''aii 



r?"^(a''/i^-5*/i^)-25°/i^^ + 5*/i°^ + 7?^^9°/i^ +0{h^). (65) 
For the khronon terms, we find that at quadratic order in the unitary gauge 

= Mi [k(-) (f + Vapf^,) - K^'fl,] . (66) 

From this expression it is clear that the notion of energy (j6ip is not well defined for spacetimes 
with radiation at infinity satisfying conditions (j38p . This is an unphysical divergence, which is 
regularized for a fiux of energy of finite duration |36| . For our purposes, it is enough to notice that 
the time variation (|64p (and hence the fiux) is well defined for these boundary conditions. Also, only 
the part of the integral quadratic in perturbations does not vanish, which implies that the previous 
expressions are enough to compute the fiux of energy at infinity. The steps to go from the previous 
formula to the final result (j40p are explained in Sec. [6j 
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